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ABSTRACT 

The conditions of encounter between interplanetary 
dust particles and the planets are investigated 
for the range of particle size where the force of 
solar radiation is not negligible in comparison 
with the force of solar gravity. The action of 
solar radiation pressure in reducing the heliocentric 
orbital speed of a small dust particle makes 
possible very low planetocentric speeds at encounter 
for an entire class of heliocentric orbits of non- 
zero eccentricity and fairly low inclination. The 
speed at encounter and the apparent radiant depend 
on the size of the dust particle for a given 
heliocentric orbit throughout the range of particle 
size being investigated by means of dust particle 
sensors mounted on spacecraft. 

INTRODUCTION 

Dust particle sensors to be flown on spacecraft in order 
to directly measure the speeds and directions of motion of small 
dust particles having masses in : ^10 gm have been developed 

during the past few years. These sensors were designed to be 
used in studying interplanetary dust particles for which the 
force of solar radiation is not negligible in comparison with 
the force of solar gravity. The range of sensitivity of the 
sensors extends into the regime of particle size where dust 
particles can be blown out of the solar system by sunlight. 
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The effects of solar radiation pressure must therefore be 
included when the measured speeds and radiants of small dust 
particles are used in computing the orbits. Evaluation of these 
effects requires investigating the conditions of encounter 
between a planet and interplanetary dust particles for which 
radiation pressure modifies the heliocentric orbital speeds. 
Conditions of encounter with radiation pressure included also 
have an important bearing on the problem of the enhanced flux 
of small dust particles measured in the vicinity of the earth 
and possibly existing also near the other planets. 

THE EARTH’S DUST BELT 

The flux of small dust particles measured in the vicinity 
of the earth through the use of rockets and satellites is higher 
than would be expected if radiation pressure is neglected and 
the orbits of small interplanetary dust particles are assumed 
to be similar to those of larger, meteor- producing meteoroids. 

This enhancement of the flux of small dust particles has been 
called a dust blanket by Beard (1959), a dust belt by Whipple 
(1961), and a dust shell by Singer (1961). The data that are 
presently available do not permit distinguishing among the various 
models envisioned by these and other investigators. A knowledge 
of the geocentric trajectories is critical in discussions of 
the nature and origin of the earth’s dust belt. Referring to 
the phenomenon as an enhancement of the flux of small dust 
particles in the vicinity of the earth carries fewer connotations 
about the presently unknown geocentric trajectories of the dust 
particles and, more importantly, permits closer contact with 
the observational data that are now available. 

The rate at which information can be obtained with a 

* 

satellite-borne dust particle sensor is increased by the enhancement 
of the flux of small dust particles near the earth. At the same 
time, the complicated geocentric trajectories responsible for 
the enhancement of the flux make the experimental determination 
of the orbits of the small dust particles more 
for meteoroids. 




The geocentric enhancement of the flux of small dust 
particles is shown in Figure 1, which is a cumulative flux 
distribution expressed in terms of the particle mass. The 
numerical data used in constructing Figure 1 are given in 
Table 1. The selected data are considered here to be representative 
of the available data, although Figure 1 is a corrected version 
of a rather outdated figure that was presented in an earlier 
paper. The distribution curves shown in Figure 1 come from 
studies of the frequency of meteorite falls, observations of 
meteors, analyses of results from photometric studies of the 
zodiacal light, and direct measurements made with dust particle 
sensors mounted on spacecraft. 

The geocentric enhancement of the flux of small dust 

particles is revealed by plotting in Figure 1 the cumulative 

flux distributions for two different regions of space. The 

data from studies of the Zodiacal light, meteors, and the 

frequency of meteorite falls are used to represent the flux of 

dust particles and meteoroids at earth’s distance from the sun 

but in regions removed from the earth. The direct measurements 

made in the vicinity of the earth are used to represent the flux 

of small dust particles near the earth. The flux of meteoroids 

is essentially the same for both regions of space, since the 

known distributions of orbits and the high speeds of encounter 

(20 to 40 km/sec, on the average, as determined by radar and 

photographic studies of meteors) preclude any appreciable geocentric 
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enhancement of the flux of meteoroids. The flux of 10 

2 -6 
particles/m /sec/2TT ster for meteoroids having masses M ~ 1.3 x 10 

observed by Elford, et al. (1964) through the use of radar 

indicates that no measurable geocentric enhancement of the flux 

occurs for the faint radar meteors. 

Comparison of the flux of small dust particles in inter- 
planetary space with the higher flux measured in the vicinity 
of the earth leads to the conclusion that a geocentric enhancement 
of the flux of small dust particles exists. The degree of 
enhancement of the flux and the particle mass at which the flux 



begins to be enhanced both depend on which of the various size 
distributions derived from studies of the zodiacal light and 
solar F-corona is used to represent the flux of small dust 
particles in interplanetary space. The size distribution derived 
by Ingham (1961) and shown in Figure 1 represents the best fit 
to the distribution of surface brightness of the zodiacal light 
that was observed by Blackwell and Ingham. This size distribution, 
if extended to smaller sizes of dust particles, is in good 
agreement with the size distribution of small dielectric spheres 
(water ice) used by Giese (1963) in calculating models for the 
zodiacal light. The photoelectric observations of the surface 
brightness and polarization of the zodiacal light reported by 
Weinberg (1964) are fitted more closely by the size distribution 
given by Giese than by any of the other size distributions that 
have been suggested. The size distribution given by Elsasser (1954) 
has subsequently been altered and was included in Figure 1 bnly to 
indicate the desirability of a smooth transition between the 
cumulative flux distribution for the zodiacal dust particles and 
the one for the radar meteors. 

The geocentric enhancement of the flux of small dust 

particles, as portrayed in Figure 1, begins to appear at a 
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particle mass between 10 and 10“ gm and reaches a maximum 
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(—10 enhancement) for dust particles having masses in the 
neighborhood of ml ~ 10” 11 gm. This dependence suggests that 
some non-gravitational force, such as that of solar radiation 
pressure, may play an important role in producing the geocentric 
enhancement of the flux of small dust particles. If so, similar 
enhancements of flux can be expected to exist in the vicinity 
of other planets. Investigation of this possibility again 
requires that radiation pressure be includedin the conditions of 
encounter between small interplanetary dust particles and the 
planets. 

The following sections are devoted to an approximate 
treatment of the conditions of encounter between the planets and 
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small interplanetary dust particles. Solar radiation pressure 
is included in the treatment, but other non-gravitational forces 
such as cospuscular drag, Coulomb drag, and Lorentz forces are 
neglected in this first order approximation. Any force which 
effectively reduces the central gravitational force exerted by 
the sun can produce effects on the motion of a small dust 
particle similar to the effect of radiation pressure. 


EQUATIONS OF MOTION AND THE POYNTING- ROBERTSON EFFECT 

Some of the equations to be used here are sufficiently well- 
known that they really need not be presented, but they do provide 
a convenient means for introducing parameters needed in the 
discussion. The scalar equations of motion given by Robertson 
(1937) or by Wyatt and Whipple (1950) form a convenient starting 
point. 


Expressed in polar coordinates (r, 0), the equations for 
the heliocentric motion of a dust particle of mass m moving in 
the combined gravitational and radiation fields of the sun of 
mass » m are 


and 
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In the foregoing, G is the constant of gravity, L© is the 
luminosity of the sun, c is the speed of light, 6 is the mass 
density of the dust particlfe, and s is the particle radius (all 
measured in cgs units). The parameter Q is the radiation 
pressure efficiency factor defined by 


'pr 


Q TTs' 
pr 


where C is the cross-sectional area of the dust particle 

pr 

effective in intercepting sunlight to produce radiation pressure 
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on the dust particle. 

The magnitude of the force of radiation is 
F r = uCr) Cp r 

where u (r) is the energy density of sunlight at a heliocentric . 
distance r. The magnitude of the force of gravity is 
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r 2 

A dimensionless parameter (3 defined as 
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will be used extensively throughout the following discussions. 

The parameter f3 is simply the magnitude of the force of radiation 
measured in units of the force of gravity. Both the flux of 
sunlight and the force of solar gravity vary inversely as the 
square of the distance from the sun, so the parameter (3 is 
independent of the heliocentric distance. Other non- gravitational 
forces could be introduced at this point by using an additional 
parameter p* which would not (in general) be independent of 
the heliocentric distance. The expression given earlier for p 
becomes, upon substitution of the appropriate numerical values 
for the various constants. 


p = 5.78 x 


io - 5 ,i- m 2 

cm 
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The equations of motion may be written also in terms of the 
dimensionless parameter (3 as 


and 
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since p and a are simply related by the expression 



Radiation pressure produces several effects on the motion 

of a small dust particle. A dust particle moving in a heliocentric 

orbit receives sunlight from a source that is moving relative 

to the dust particle, so the incident sunlight undergoes 

aberration through a small angle 

r d9 
* ~ c dt 

which can be neglected in a first approximation. The component • 
of motion of the dust particle along the position vector of the 
dust particle causes sunlight of frequency v to be shifted 
through the Doppler effect to a new frequency v* given by 

v' = v(l - i 

The radial motion of the dust particle and the shift in frequency 
of the incident sunlight produce a corresponding change in the 
energy density of the intercepted radiation and, hence, the 
radiation pressure exerted on the dust particle. The terms 
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represent the radial and tangential components, respectively, 
of the acceleration produced by this drag force. The equations 
of motion given previously are correct to the first order in 
the ratio of the speed of the dust particle to the speed of light. 
The problem was first formulated correctly by Robertson (1937) 
who used a relativistic treatment in order to derive the 
generalized equations of motion. The action of sunlight in 
effectively reducing the central gravitational force of attraction 
by the sun and in giving rise to the non-conservative drag 
force is called the Poynting- Robertson effect. 

The energy and angular' momentum no longer constitute 
constants of motion because of the presence of the non-conservative 
drag force introduced by radiation pressure* Integration of 
the angular equation of motion yields 


h = r^ cM3 = ho - jig 9 
dt c 

where h is the angular momentum (per unit mass) of the dust 
particle and ho is the value of h at some epoch, t =* t Q . As 
time t (or 0) increases, the angular momentum decreases. The 
dust particle loses angular momentum and energy of motion, 
causing it to spiral toward the sun. The dynamical evolution 
of the orbits of meteoroids under the Poynting-Robertson effect 
has been considered by Wyatt and Whipple (1950) for eBLiptic 
orbits as well as for the circular orbits considered by Robertson 
(1937). 

Temporarily dropping the small terms in 
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from the equations of motion leaves 
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The radial equation of motion now displays only the action of 
solar radiation pressure in effectively reducing the central 
force of gravity. A value 3 > 1 makes the right-hand side of 
the radial equation of motion positive and causes the central 
force to become one of repulsion. A dust particle for which 
(3 > 1 will be blown out of the solar system by radiation pressure, 
while one for which p = 1 will not be subject to any acceleration 
by the sun. 

Integration of these approximate equations of motion yields 
the orbital speed v of the dust particle as 

v 2 - ua-fs) cf - b 

where a is the semi-major axis of the orbit of the dust particle. 
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Comparison with the analogous equation 



for a dust particle which is sufficiently large that p = 0 shows 
that radiation pressure reduces the orbital speed of a small 
dust particle below that for the purely gravitational case. A 
small dust particle moving in a planet-like orbit moves more 
slowly than the planet and therefore has a non-zero speed 
relative to the planet. This effect was noted by Poynting in 
1912, soon after the existence of radiation pressure had been 
experimentally verified. 

An assumption commonly used in simplified treatments of 
this type is that a dust particle is perfectly absorbing and 
absorbs sunlight over an area equal to the geometric cross- 
sectional area of the particle. Values of (3 computed under this 
assumption (Q pr = 1) are shown graphically in Figure 2 for 
various values of the mass density 6 and particle radius s. 

As was noted previously, a value 3 >1 means that the dust 
particle will be blown out of the solar system by radiation 
pressure. The values p = 1 occur in Figure 2 for particle radii 
s ~ In which are comparable to the wavelength of sunlight 
(~0.55 m.). The assumption that classical optics could be used 
to compute P fails. The radiation pressure efficiency factor 
Q pi ; (and, hence, p) is one of the parameters which can be 
computed through use of the Mie theory of light-scattering by 
spherical particles of arbitrary size and known index of refraction. 
It is unlikely that interplanetary dust particles are isotropic 
homogeneous spheres, but one presently has little choice in 
computing the values of p other than to assume a composition 
for which the index of refraction is known and to use the Mie 

t 

theory for computing Q pr . 

Further treatment of the rather complicated relations among 
particle mass or size, mass density, composition, structure, 
shape, and the parameter p are beyond the scope of this work. 

An extensive treatment of the Mie theory can be found in the 
book by van de Hulst (1957) . The important points in the present 



discussion can be adequately illustrated by taking Q = 1 and 
using Figure 2 for relating (3 to the physical parameters of the 
dust particles, A scale of particle mass computed for a mass 
density of 6 ==2,5 gm/cm has been added along the top border 
in Figure 2 for this purpose. 

The direct measurements apply for dust particles having 
-7 

masses m ^ 10 gm. Reference to Figure 2 shows that such dust 
particles have values of f3 ^ 0.1, approximately, for any of the 
values of mass density 6 that are represented. A value (3 = 0.01 
might be considered to be negligible in comparison with unity, 
but such a value applies only for the upper end of the range 
of particle mass encompassed by the direct measurements. The 
majority of the direct measurements apply for dust particles 
having masses between about 10” gm and 10 gm. A dust 
particle having a mass m = 10” ^ gm and a mass density 6 =2.5 

3 

gm/cm has a radius of about 1 p. and, from Figure 2, a value of 
(3 a 0.23, which is hardly negligible compared to unity. 

The dependence of the speed at encounter (entry of the dust 
particle into the sphere of influence of a planet) on the value 

_7 

of (3 as the mass of the dust particle is varied from m ~ 10 
gm down to any limits ((3 = 1) set by radiation pressure is to be 
emphasized. The treatment is approximate because of the simplifying 
assumptions made in relating |3 to the particle mass or the 
particle size. It is also approximate to the extent that the 
heliocentric motion of a dust particle is followed by using 
the two-body equations of motion presented earlier with radiation 
pressure included. As the dust particle enters the planet*s 
sphere of influence, gravitational control is transferred to 
the planet, and the motion of the dust particle in the region of 
the planet is treated as a two-body problem. The small tangential 
drag force produced by radiation pressure can be neglected during 
an encounter but should be included in the motion of a dust 
particle between encounters. 
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The orbital speed of a dust particle at encounter is given 
by 

v 2 = n(l-p) (2 - i) 

in which the semi-major axis a of the heliocentric orbit of the 
dust particle is measured in units of the semi-major axis a^ of 
the planet’s orbit (assumed to be circular) as 



The notation to be used here is that employed by Opik (1951) . 

The relative speed can then be found by vector addition of the 
velocity of the dust particle and the velocity of the planet. 
Alternatively, and much more simply, the problem can be formulated 
as a Jacobi three-body problem with radiation pressure included. 
The speed U of the dust particle at encounter (measured relative 
to the planet and expressed in units of the circular speed v c of 
the planet) is given by 

U 2 = (J-) 2 =1-2 -/l-V" -W 1-e 2 ) cos i + (1 - p)(2 - b 
v c A 


in which e is the eccentricity of the heliocentric orbit of the 
dust particle and i is the inclination of the heliocentric orbit 
of the dust particle measured relative to the plane of the 
planet’s orbit. Setting i = 0° and replacing A by Q where 

Q = A(l-e) , q = a(l-e), and Q = ■§ 

a o 


yields 


U 2 = 


1-2 


1-P IQ \l+e) + (1-p) (2 - - ( --^ ) . 


The foregoing expression reduces to 

U 2 = 1 - 2 ~f~ L+e ~ + (1-p) (1+e) 


for encounter at the perihelion passage (Q = 1) of the dust 
particle. This expression is quadratic and has the root 


U = 1 
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Now, U = 0 for (3 = P Q where P Q and e are related by the expression 

Po = 1 — 1 . 

TFe 

The ranges of p Q and e for which U = 0 can occur when radiation 
pressure is included are 

0 £ P Q . < 0.5, 0 £ e < 1, (elliptic orbits) 

Po = 0.5, e ** 1, (parabolic orbits) 

0.5 < P Q < 1, e>l, (hyperbolic orbits). 

Circular heliocentric orbits (e = 0) give U = 0 only for p = 0 

which applies to dust particles considerably larger than those 

studied through the direct measurements technique. Figure 3 

shows U (or v/v.) as a function of p for circular orbits and 

c 

for several values of the inclination. The ordinate on the 

right will be explained and used later. It should be noted 

that small dust particles having p ~ 0.01 and planet- like orbits 

do not have particularly low speeds of encounter (even for zero 

inclination) when the effect of radiation pressure on the orbital 

speed of a dust particle is included. 

Figures 4a and 4b show U (or v/v c ) as a function of P 

for two values of eccentricity e arbitrarily chosen such that 

P Q falls in the range of particle size (P ~0.01) encompassed by 

the direct measurements. Figures 3, 4a, and 4b show that very 

low speeds of encounter between planets and small dust particles 

can occur for dust particles having non-circular heliocentric 

orbits of fairly low inclination. Further calculations show 

that the singular solutions shown in Figures 4a and 4b are but 

two of an entire class of orbits for which small dust particles 

can have low speeds at encounter. Figure 5 shows the (A - "*", e) 

# 

domain in which the clear areas contain values of A and e for 
which planetary encounters are possible. Representative values 
of p 0 are noted along the line for perihelion intercepts. The 
class of singular solutions which give U * 0 for i = 0° forms 
a line singularity. An eccentricity e £ 0.01 is required in 

order for P Q to fall in the range (p * 0.01, m^ 10-7 gm) 



encompassed by the direct measurements. The entire class of 
orbits for which i ^ 5 °, e ^ 0.01 provides small dust particles 
which have low speeds at encounter provided the encounter occurs 
near the perihelion passage of a dust particle. The action of 
radiation pressure in reducing the perihelion distance of a dust 
particle through the Poynting-Robertson effect further serves 
to continuously move orbits of low inclination into position 
for the low speed encounters to occur. The line singularity 
appearing at Q = 1 in the (A“\ e) domain for i = 0° presents 
no problem regarding an "inf ini te" enhancement of the flux in 
the vicinity of a planet, because the number of dust particles 
having such orbits if finite. The orbits of primary interest 
in the problem of the enhanced flux near a planet are those in 
the neighborhood of the line Q = 1, 0,1 ^ e < 1, i » 0°, in 
the (A" 1 , e) domain. 

The component of U normal to the plane of the planetary orbit 
is given by 

U z = i “fT-P ' - /A(l-e 2 ) sin i, 

so increasing i from 0° adds a non-vanishing component to U. 
fairly rapidly through the dependence of U ^ on sin i. The 
factor broadens the range of i over which a given small 

value of Ug can be obtained. The component of U along the 
heliocentric position vector of a dust particle is 

U x - ± '/IT £(2 - i) - A(l-e 2 )j 1/2 . 

A non-vanishing component is added to U fairly rapidly as Q is 

reduced below 1. Again, the factor -fl - (3 reduces the magnitude 

of the non-vanishing component. Evaluation of the degree of 

enhancement of the flux of small dust particles which will occur 
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near a planet therefore requires an integration over the orbital 
elements a, e, and i as well as over a size distribution or a 
distribution of the values of p for the dust particles encountered 
by the planet. The writer has not yet been successful in 
performing this integration. 
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The action of radiation pressure in reducing the orbital 
speed of a small dust particle also affects the apparent 
radiant of the dust particle. Varying the size of a dust particle 
so that the value of p crosses a singularity of the type shown 
for i = 0° in Figures 4a and 4b causes the radiant of the dust 
particle to shift by 180°. Dust particles having orbits with 
e > 0, i => 0°, and P < P Q appear to overtake the planet, while 
those in similar orbits but having P > P Q appear to be overtaken 
by the planet. If a planet encounters a cloud of dust particles 
all having similar orbits of low inclination and a distribution 
of values of p encompassing the value of p„ appropriate for 
the eccentricity of the orbits, the dust particles can appear 
to enter the planet’s sphere of influence from all directions. 

Such a cloud of dust particles would, of course, be in the 
process of dispersing because of the differential effects of 
radiation pressure on the heliocentric orbits of the dust 
particles. 

The effects of radiation pressure on the conditions of 
encounter that were described in the preceding paragraphs 
apply also to the small dust particles detected by a dust 
particle sensor mounted on an interplanetary probe such as 
Mariner IV. The effects of radiation pressure on the speeds 
and radiants of the dust particles detected by such an instrument 
need to be investigated in greater detail. 

The effect of radiation pressure on the accretion process 
has some implications (the importance of which has not yet been 
investigated) regarding the formation of planetesimals in a 
rotating circumstellar nebula out of which a star forms. As a 
newly-formed star becomes luminous and clears away the gas 
in the nebula, the small du&t particles which are not blown 
away by radiation pressure assume Keplerian orbits about the 
star. The accretion process may then change because of the 
differential effects of radiation pressure on the planetesimals 
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and the smaller dust particles being accreted by the planet- 
esimals. The solar system is in a stage somewhat analogous 
to the latter. The nebular gas is gone, so the planets are now 
accreting small dust particles which have their heliocentric 
motion modified by radiation pressure and the radial flow of 
the solar wind. 

GRAVITATIONAL ENHANCEMENT OF THE FLUX OF INTERPLANETARY DUST PARTICLES 

The geocentric enhancement of the flux of small dust 
particles has been the subject of considerable discussion both 
before and after its existence was recognized in the results 
obtained with dust particle sensors flown on rockets and satellites. 
Whipple (1961) suggested that the dust particles were the product 
of meteoroidal impacts on the moon. The majority of the remaining 
proposed explanations for the existence of the enhanced flux 
basically depend on the gravitational concentration of interplanetary 
dust particles streaming at low speed into the vicinity of the 
earth. No attempt is made here to review the numerous variations 
on this theme, because the effects of radiation pressure on the 
conditions of encounter have been neglected. 

The gravitational enhancement of the flux of interplanetary 
dust particles entering the vicinity of the earth has been 
treated through the use of Liouville’s theorem by Singer (1956, 

1961) and by Shelton, et al. (1964). The three-body problem 
(consisting of the sun, the earth, and a dust particle) was 
treated as a succession of two-body approximations. Singer (1961) 
assumed that the velocities of the dust particles were isotropic 
at infinity (entry into the earth’s gravitational field). 

The degree of enhancement of the flux depends on the average 
speed at which dust particles enter the vicinity of a planet. 

The requirements on the speed of encounter can be illustrated 
by a very brief treatment using Kiouville’s theorem and following 
closely the expositions given by Singer (1961) and by Shelton, 
et al. (1964). The general treatment applies for any planet, 
with the degree of enhancement of the flux depending rather 
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sensitivity on the distributions of a, i, and p for the dust 

particles encountered by the planet. The remainder of the 

present discussion is limited, however, to the problem of the 

geocentric enhancement of the flux of dust particles. 

The densities at two points (r^, o and (r*, v*) on a 

g g g g 

dynamic trajectory in a six-dimensional phase space composed 
of the geocentric position and velocity vectors r and v are 
related through Liouville f s theorem by 
n <^, T g ) = n (r>, v*) 

if no particles are lost. The density of dust particles in real 
space at a geocentric distance r^ is obtained by integrating 
n(r , v ) over velocity space to get 

D D . w 

p (r g )=( n(r g , ^)dv g 

J g 

The scalar flux of dust particles at a geocentric distance r g 
is 


M'g) " P(r g ) v g (r g ) 


taking r», v», 
the earth’s gravitational field leads to 

■ f [i + v i ( V ]1/2 
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and to apply when the dust particles enter 


and 


<V 




< V 


= f 




[i + v l (r g ) ] 
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as the density enhancement and flux enhancement, respectively, 
at a geocentric distance r . The factor f allows for those dust 

o 

particles which are absent from the receding flux because they 
collided with the earth. The algebraic expression for the 
factor f is given by Shelton, et al. (1954). The numerical 
values of f range between 0.5 and 1. The speed at infinity 
v corresponds to the speed at encounter v used earlier in the 
discussion of the conditions of encounter. 
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The conservation of energy of geocentric motion yields 

v^(r) = v‘ 
g g 


v e<V 


where 


v (r 
e ' g 






g 


is the speed of escape from the earth at a distance r 


g 


At the 


altitudes of the near-earth satellites (r^ ~ 7000 km), the speed 
of escape is 


v ~ 10 km/sec. 
e 

For speeds at encounter 


« v ~ 10 km/sec, 
e 


one has 

A 


( V ~ v e (r g ) 


so that the density and flux enhancements reduce to 

(r g } - f v e 


g' 


Oo 


and 

< < ' r g ) ~ 1 T e (r g ) . 

~T~T 2 

4 

An enhancement of the flux by a factor ~ 10 therefore 
requires ~ 0.1 km/sec. The concept of the sphere of 
influence inherent in the simplified treatment using Liouville*s 
theorem encounters trouble, for 0.1 km/sec is not small compared 
to the speed of escape at the distance of the boundary of the sphere 
of influence for the earth. A realistic treatment of the problem 
therefore requires numerically integrating the three-body equations 
of motion with radiation pressure included. 

The orbits for which very low speeds at encounter are possible 
for dust particles of various sizes were discussed earlier and 
represented in Figures 3, 4a, 4b, and 5. The speed at encounter 
v measured in units of the speed of escape v 0 at the altitudes 
of the near-earth satellites is given as the ordinate on the 
right in Figures 3, 4a, and 4b. 
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Singer (1961) assumed that relatively large dust particles 

(s ~ 200ia) have almost circular heliocentric orbits of very low 

inclination (e ^0, i ^0°) in order to justify taking v _ =* 1 km/sec, 

2 ^ 

which leads to a flux enhancement ~10 . A subclass of these 

orbits was investigated by Dole (1962), who numerically integrated 

the geocentric trajectories of dust particles having initially 

circular heliocentric orbits of zero inclination (e = 0, i = 0°). 

Southworth (1963) has criticized the two-dimensional model used 

by Dole and has argued that the fraction of zodiacal dust particles 

having such heliocentric orbits must be quite small. 

A dust particle having a radius s ~ 200(jl and a mass density 
3 -3 

of 1 gm/cm has a value of P~ 2.4 x 10 , which justifies the 

neglect of radiation pressure by Singer (1961). But the same 

-5 

dust particle has a mass m ~ 3.3 x 10 gm, which removes it from 

the regime of particle mass covered by the direct measurements 

and places it in the regime of radar meteors. The available direct 

< -7 

measurements apply for dust particles having masses m ~ 10 gm, 
with the maximum enhancement of the flux occurring for dust 
particles having masses in the neighborhood of m ~ 10” ^ gm. A 
dust particle having a mass m = 10 -11 gm and a mass density r6 = 2.5 
gm/cm has a radius s ~ 1 |i. The value of p is about 0.23, 
according to Figure 2. The speed of encounter for such a dust 
particle having an earth-like orbit is about 3.6 km/sec. Taking 
this as an average speed yields a flux enhancement of only 8, 

4 

which is very much less than the factor ~ 10 given by the direct 
measurements. 

The shortcomings of the treatments by Singer (1961) and Dole 
(1962) are shared by numerous other proposed explanations for 
the geocentric enhancement of the flux of small dust particles 
in which radiation pressure Was neglected. The work of Dole 
was mentioned here because such numerical integrations should be 
extended to a three-dimensional model and to the case in which 
radiation pressure is not neglected. The formulation of the 
problem of the gravitational enhancement of flux in terms Liouville*s 
theorem that was given by Singer (1961) is useful, provided two 
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important changes are made. These are (1) radiation pressure 
should be included so that the discussion can be directed to the 
range of particle size for which the measured geocentric enhancement 
of the flux occurs, and (2) the assumption that the small dust 
particles which produce the zodiacal light have predominantly 
circular heliocentric orbits of very low inclination should be 
avoided. The assumption that small dust particles having orbits 
of fairly low inclination contribute most to the enhancement of 
the flux near a planet cannot be avoided, but the restriction on 
the range of the inclinatio n can be relaxed slightly because of 
the presence of the term ~fl~ - 0 in the expression for U z . The 
low speed encounters occur only for non-circular orbits of fairly 
low inclination when dust particles of the sizes involved in the 
measured enhancement of the flux are considered. 

The remaining unsolved problem consists of integrating over 
the unknown distribution of orbits and the rather poorly determined 
size distribution for the small interplanetary dust particles in 
an attempt to ascertain whether appreciable planetocentric 
enhancements of the flux can occur directly. The boundary 
conditions which must be met by any distributions of orbits and 
size distributions assumed for use in such an integration include 
the observed distribution of surface brightness and polarization 
of the zodiacal light over the celestial sphere. 

SUMMARY 

The conditions of encounter between planets and small 
interplanetary dust particles have been investigated in an approximate 
manner,. Solar radiation pressure has been included in the conditions 
of encounter, and the dependence of the speed at encounter on 
the size of a dust particle' has been emphasized. 

The low speeds of encounter which occur near perihelion 
passage encompass an . entire class of non-circular heliocentric 
orbits of fairly low inclination when dust particles of the sizes 
being studied through the use of spacecraft are considered. This 
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class of orbits provides dust particles which can encounter the 
planets at the very low speeds necessary in order for a gravitational 
enhancement of the flux Id occur. Evaluation of the effectiveness 
of this class of orbits in directly producing the enhanced flux 
of small dust particles measured in the vicinity of the earth is 
made difficult by the present lack of knowledge about the orbits 
of small dust particles both in the vicinity of the earth and in 
interplanetary space. The source of the small dust particles 
and the mechanism for the enhancement of the flux measured in 
the vicinity of the earth through the use of spacecraft still 
remain to be satisfactorily explained. 


Figure 


Figure 


Figure 


Figure 


Figure 


CAPTIONS FOR FIGURES 


1. Cumulative flux distribution compiled from the data 
in Table I to show the enhancement of the flux of 
small dust particles in the vicinity of the earth. 

2. Dependence of the ratio ( 3 ) of the force of radiation 
to the force of gravity on the size and mass density 
of a dust particle. 

3. Dependence on particle size (or 3 ) of the speed at 
encounter for dust particles and a planet when both 
are moving in circular heliocentric orbits. 

4a, 4b. Two examples of the dependence on particle 

size (or 3 ) of the speed at encounter for a planet 
moving in a circular orbit and dust particles 
moving in elliptic heliocentric orbits. 

5. The (A \ e) domain showing the orbits for which 
planetary encounters are possible and the line 
singularity next to which are clustered the orbits 
that yield very low speeds of encounter for small 
dust particles in orbits of low inclination. 
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